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(Dated: August 19, 2005) 

We present the results ol a search for the effects of large extra spatial dimensions in pp collisions 
at ^/s — 1.96 TeV in events containing a pair of energetic muons. The data correspond to 246 
pb _1 of integrated luminosity collected by the D0 experiment at the Fermilab Tevatron Collider. 
Good agreement with the expected background was found, yielding no evidence for large extra 
dimensions. We set 95% C.L. lower limits on the fundamental Planck scale between 0.85 TeV and 
1.27 TeV within several formalisms. These are the most stringent limits achieved in the dimuon 
channel to date. 
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In their 1998 paper, Arkani-Hamed, Dimopoulos, and 
Dvali (ADD) suggested that the seemingly unreachable 
Planck energy scale (conventionally thought to be Mpi ~ 
10 19 GeV) may be in fact much lower, i.e. within the 
reach of current and planned future colliders 0. They 
postulated that the standard model (SM) particles and 
gauge interactions are confined to a three-dimensional 
"brane" embedded in a "multiverse" , which consists of 
the three standard plus n additional compact spatial 
dimensions. However, gravitons in this framework can 
propagate in the entire multiverse. The gravitons propa- 
gating in compact extra dimensions appear as a tower of 
Kaluza-Klein (KK) excited modes from the point of view 
of the SM brane. Furthermore, the radius of compactifi- 
cation (R) of extra dimensions in the ADD model is much 
larger than either the Planck or electroweak length, and 
may be as large as ~ 1 mm. Since gravitons are free 
to propagate in these large extra dimensions, the gravi- 
tational interaction would appear suppressed on the SM 
brane, due to the extra volume gravity permeates. Con- 
sequently, while the apparent Planck scale is ~ 10 19 GeV, 
with respect to the 3 + n-dimensional space, the funda- 
mental Planck scale (Ms) can be as low as ~ 1 TeV, thus 
eliminating the hierarchy problem of the SM. 

The phenomenological consequences of the ADD model 
have been a subject of intense study in recent years. For 
a review of the possible effects of large extra dimensions, 
ranging from modification of Newtonian gravity to black 
hole production at future colliders, see, e.g., Ref. 

In this Letter, we describe a search for the effects of 
large extra dimensions via virtual Kaluza-Klein graviton 
(Gkk) exchange in pp collisions resulting in the dimuon 
final state. Technically, virtual graviton effects are sen- 
sitive to the ultraviolet cutoff required to keep the di- 
vergent sum over the KK states finite 0, 0, |g, rather 
than the fundamental Planck scale. As the two scales 
are expected to be closely related, we do not distinguish 
between them in this analysis. The search is based on 
246 ± 16 pb" 1 of data collected in 2002-2004 by the D0 
detector operating at the Fermilab Tevatron Collider at 
y/s = 1.96 TeV. We used the method of Refs. 0, 0, in 
which the dilepton invariant mass (M) and the cosine of 
the scattering angle (cos 8*) in the dilepton center of mass 
frame are analyzed simultaneously for the effects of large 
extra dimensions. This is the first search for large extra 
dimensions in the dimuon channel at a hadron collider. 
Previous searches for virtual graviton effects in various 
channels at LEP, HERA, and the Tevatron are reviewed 
in detail in Ref. [3j. 

The dimuon production cross_section in the presence 
of extra dimensions is given by 



dMdcos t 



= fsM + /int VG + /KK Vg> 



(1) 



effects of large extra dimensions are parametrized via a 
single variable rjc = F/Mg, where T is a dimensionless 
parameter of order unity. Three different formalisms for 
T are explored in this analysis: 



T = 1 

[Giudice-Rattazzi- Wells (GRW) 0]; 



[Han-Lykken-Zhang (HLZ) 0]; 



n = 2 
n > 2 



(2) 



(3) 



2A 2 
T = — = ±- 

7T 7T 



(4) 



where fsM, /int, and /kk are functions of (M, cosd*) 
and denote the SM, interference, and Gkk terms. The 



(Hewett 0). In Eq. J3J A = ±1 indicates whether vir- 
tual graviton exchange interferes constructively or de- 
structively with SM processes. While virtual graviton 
exchange does not depend strongly on n (the number of 
extra dimensions), the HLZ formalism for T does explore 
this dependence explicitly. 

The D0 detector and its data acquisition system are 
described in detail elsewhere [j| . Here we give a brief de- 
scription of the components used in the analysis. At the 
center of the D0 detector is the central-tracking system, 
which consists of a silicon microstrip tracker and a central 
fiber tracker. Both trackers are located within a 2 T ax- 
ial magnetic field |9j , with designs optimized for tracking 
and vertexing at pseudorapidities \r]d\ < 3. The pseudo- 
rapidity rj is defined as — In (tan |) , where 9 is the polar 
angle w.r.t. the proton-beam direction, as measured from 
the interaction vertex. We also define which is the 
pscudorapidity based on the polar angle measured from 
the geometric center of the detector. The nearly hermetic 
uranium/liquid-argon calorimeter is used to measure en- 
ergies of electrons, photons, and hadrons The muon 
system covers |?7d| < 2 and consists of a layer of track- 
ing detectors and scintillation trigger counters in front of 
1.8 T iron toroidal magnets, followed by two more similar 
layers of detectors outside the toroids 01 • Luminosity 
is measured using plastic scintillator arrays placed in the 
large \rj d \ (forward) regions of the detector. 

The sample of candidate events used in the search was 
collected with a set of triggers that require either one 
or two muon candidates in the muon system. After de- 
tailed event reconstruction an event must contain at least 
two muon candidates, each matched with a track in the 
central tracker, which is used for muon momentum mea- 
surement. The matching tracks were required to have 
transverse momenta pt > 15 GeV, at least one hit in the 
silicon microstrip tracker, and at least nine hits in the 
central fiber tracker. The latter two requirements ensure 
reliable momentum measurement, especially at high px- 
To reduce background from cosmic rays, we introduced 
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additional criteria. Since cosmic muons are not corre- 
lated with the beam crossing, we required the muon ar- 
rival time, as measured in the muon system scintillation 
counters, to be within 10 ns (four standard deviations) of 
the expected arrival time for an highly-rclativistic parti- 
cle produced in a pp collision in the center of the detec- 
tor Furthermore, dimuon events that originate from 
a cosmic muon are back-to-back in rj. This is because the 
same cosmic muon is reconstructed twice in the event, 
once when entering the detector and once when leaving 
it. Consequently, dimuon events from cosmic rays should 
have rji + 7/2 ~ 0. In true dimuon events originating from 
pp collisions, the two muons are generally not back-to- 
back in 77 due to a longitudinal boost of the dimuon sys- 
tem. Therefore the sum 7/1 + 772 was required to be away 
from zero by at least five standard deviations of the Gaus- 
sian distribution as observed in a cosmic ray sample. The 
signal efficiency for this selection is (99 ± 1)%. After the 
above selections the cosmic muon contamination in the 
candidate sample is negligible. 

Muons from graviton decay are expected to be iso- 
lated from other energetic particles or jets. In contrast, 
high pt muons from b and c quark decays tend to be 
non-isolated. To reduce this copious background, we re- 
quired that each muon have: (i) J2tz=q s(Pt) < 2.5 GeV, 
where X)k=o s(pt) is the scalar sum of the px of all ad- 
ditional tracks contained within a cone of radius 1Z = 
y/ (A77) 2 + (A0) 2 = 0.5 ^2 centered on the muon track; 
and (ii) J2k=oa( e t) ~ J2tz=o.i( e t) < 2.5 GeV, where 
the Et's are the transverse energies of the calorimeter 
cells within the respective cones centered on the muon 
track. 

After all the selections are applied the efficiency per 
muon is (80 ± 4)%, as measured with Z — > /_t/x events. 
This includes the efficiency for muon and track recon- 
struction, track matching, number of tracker hits, and 
cosmic ray muon vetos, as well as the isolation selections. 

Because this analysis focuses on very high-p^ objects, 
the most problematic background is that from mismea- 
sured Drell-Yan (DY) events that appear to have very 
high mass. This occurs because if a muon's px is mis- 
measured toward a higher value it tends to be grossly 
mismeasured high, and thus the reconstructed mass also 
tends to be much higher than it actually is. This is due 
to the resolution of the central tracker, which is approxi- 
mately Gaussian in 1/pr with a typical a = 0.00272. We 
reduce the effect by scaling the px of each muon track to 
a weighted average based on the original track — mea- 
surement and its uncertainty. For example, in tne high- 
est mass event the original muon px's were 250 GeV and 
1000 GeV, while after p^-fixing both p-r's became 400 
GeV. This procedure assumes the two muons' transverse 
momenta should be equal. For high mass objects such 
as a KK graviton the equal momenta assumption would 
be accurate. The following equation illustrates the pro- 



cedure: 

J_ = J_ = \wi/p Tl +€W 2 /PT 2 \ , g s 
p' Ti p' T2 W1+W2 

where pt x 2 are the original pt's of the two muons, W\ 2 — 
1/ct 2 (1/pti i2 ) are the Gaussian weights, and a(l/pT li2 ) 
is the uncertainty on 1 /pti 2 measured from the shape 
of the Z boson peak. In the equation above e = +1 for 
muon pairs with opposite charge and e = — 1 for same 
charge muon pairs. This search does not require that 
the two muons have opposite charges, because the effi- 
ciency of such a requirement degrades quickly at high 
masses and the requirement does not reduce the already 
low background there. 

Based on the values of l/p' T and the original angular 
information, a new momentum was calculated for the 
two muons in the event. The same procedure was also 
applied to the simulated signal and background. After 
the reweighting was done we required the dimuon mass 
M to be greater than 50 GeV, which resulted in a final 
candidate sample of 17,128 events, most of them in the 
vicinity of the Z boson peak. 

We modeled the effects of large extra dimensions via 
the parton-level leading-order (LO) Monte Carlo (MC) 
generator of Ref. augmented with a parametric sim- 
ulation of the D0 detector. The simulation takes into 
account detector acceptance and resolution for muons. 
The generator includes effects of initial state radiation as 
described below, and different parton distribution func- 
tions (PDF). We used the leading order CTEQ5L 
PDF to estimate the nominal prediction. The param- 
eters of the detector model were tuned using Z — ► fifi 
events. The simulation includes SM DY contributions 
(Z/7*), Kaluza-Klein graviton exchange diagrams, and 
their interference. 

Since the MC generator contains only the LO parton- 
level processes, we modeled next-to-leading order (NLO) 
effects by adding a transverse momentum to the dimuon 
system. The model is based on the transverse momen- 
tum spectrum of dielectron candidates observed in the 
calorimeter, which provides better high energy resolution 
capability than that of the central tracking system. Since 
the parton-level cross section is calculated at LO, we ac- 
counted for the NLO enhancement in the SM background 
by scaling the cross sections by a constant if -factor of 
1.34 Q. We assigned a ±5% systematic uncertainty on 
the value of the if -factor to account for its mass depen- 
dence. We assumed the same constant if -factor for the 
extra dimensions signal. Recent NLO calculations of the 
virtual graviton exchange cross section 0] showed that 
such a choice of signal if -factor is a reasonable assump- 
tion. That study calculated the NLO if-factor to be 1.3 
for masses around 500 GeV. 

The main SM source of isolated dimuons is DY pro- 
duction, which was modeled via the MC discussed above. 
Other SM sources (such as bb, Zj, WW, Z — > tt, and 
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FIG. 1: Comparison of data (circles with error bars) and SM predictions (histogram) in pr, r/d, and 
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ti production) are negligible, as they either have small 
cross sections compared to that for DY or are eliminated 
by the event selection. The SM DY background predic- 
tion reproduces the main kinematic characteristics of the 
candidate sample, as illustrated in Fig. ^ 

Figure |3 shows the two-dimensional distribution in M 
vs. |cos 6** | for the SM background (other backgrounds 
are negligible and ignored), the sum of the background 
and an extra dimensions signal of r/c = 3 TeV -4 , and 
the data. The data agree with the background predic- 
tion and do not exhibit evidence for large extra dimen- 
sions, which would produce an excess of events at high 
mass. The two highest-mass events, while having prop- 
erties typical for the signal, are still in good agreement 
with the SM predictions alone. To further illustrate the 
agreement between the SM background and data, Fig. [3] 
shows the one-dimensional mass distribution. For refer- 
ence, the background prediction is 4 events for masses 
greater than 400 GeV and we see 3 events in the data. 

We set limits on the fundamental Planck scale Ms via 
a Bayesian fit to the data with the full signal cross section 
given by Eq. JJJ in the entire (M, |cos6>*|) plane shown 
in Fig. The fit parameter t\q was assumed to have a 
flat prior distribution. Systematic uncertainties on signal 
and background were accounted for in the fit and include 
if-factor shape (5%), the modeling of the pt smearing 
in the MC (6%), the dependence on the choice of PDF 
(5%), pt dependence of the muon efficiency (5%), and 
the MC-to-data normalization fit (1%). The latter uncer- 
tainty accounts for the fact that we used Z — > /i/i events 
in the signal sample to find the overall normalization for 
the MC predictions, which reduced significantly the over- 
all uncertainty on the product of the efficiencies and the 
integrated luminosity. This normalization was done in 
the Z boson mass peak region from 50 GeV to 120 GeV. 

The best estimates of the parameter t\g found in this 



analysis are: 

VG = 0.00 ±g gg TeV" 4 ( VG > 0) (6) 
t]g = -0.36 ±0.35 TeV~ 4 (r, G < 0), (7) 

which are fully consistent with the SM value of i]g = 0. 
From this the one-sided 95% C.L. limits on rjo are: 

TjQ < 0.76 TeV" 4 (t/g > 0) (8) 
t) G > -0.84 TeV" 4 (tig < 0). (9) 

Our results are in good agreement with the expected sen- 
sitivity, as obtained by an ensemble of MC trial exper- 
iments (0.76 TeV" 4 for i lG > 0). The use of both the 
mass and angular variables in the fit allowed for w 7% 
improvement in the sensitivity to tig- 

We express these results in terms of limits on the fun- 
damental Planck scale within the three formalisms of 
Eqs. I0H0- In the formalism of Hewett 0, both signs 
of tig are possible and therefore both limits, B and J§J , 
are relevant. In the other two formalisms [4j,|5j, t\g IS al- 
ways positive, and only the first limit is relevant. For the 
HLZ formalism, the case of n — 2 is special since J 7 , and 
therefore i]g, depends on M 2 . To relate tig to Ms for 
n = 2, we used an average M 2 for the Gkk term at the 
Tevatron of (0.64 TeV) 2 |7j. The limits are summarized 
in Table IJ 

TABLE I: Lower limits at the 95% C.L. on the fundamental 
Planck scale, Ms, in TeV. 

GRW HLZ Hewett 

n—2 n=3 n=4 n=5 n=6 n—7 X = +1 A = — 1 
1.07 1.09 1.27 1.07 0.97 0.90 0.85 0.96 0.93 



In summary, we have performed the first search for 
large extra spatial dimensions in the dimuon channel at 
hadron colliders by looking for effects of virtual Kaluza- 
Klein gravitons. We found no evidence for large extra 
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FIG. 2: Two-dimensional distributions in the dimuon mass vs. |cos#*| for: (a) SM background, (b) the sum of the SM and 
large extra dimensions contributions for t\g — 3 TeV -4 , and (c) data. 
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dimensions in this channel with « 250 pb -1 of data col- 
lected in Run II of the Fermilab Tevatron Collider. We 
set a 95% C.L. upper limit of 0.76 TeV~ 4 on the param- 
eter t]g (for t\q > 0) that describes the strength of the 
extra dimensions effects. This result corresponds to lim- 
its on the fundamental Planck scale ranging between 0.85 
and 1.27 TeV for several formalisms and numbers of large 
extra dimensions. For comparison in Run I D0 placed 
a limit on Ms, in dielectron plus diphoton production, 
of 1.1 TeV in Hewett's A = +1 formalism, while LEP's 
DELPHI experiment placed a limit on Ms in dimuon pro- 



duction of 0.73 TeV in Hewett's A = +1 formalism |3j. 
The limits from this analysis represent the most restric- 
tive achieved in the dimuon channel to date. Results 
presented here also represent the most precise test of 
high mass SM Drell-Yan production in this channel at 
a hadron collider. 
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